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Multiple Reflections and Diffuse Scattering in Bragg Scattering at Optical Lattices 
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We study Bragg scattering at ID atomic lattices. Cold atoms are confined by optical dipole 
forces at the antinodes of a standing wave generated inside a laser-driven cavity. The atoms arrange 
themselves into an array of lens-shaped layers located at the antinodes of the standing wave. Light 
incident on this array at a well-defined angle is partially Bragg-reflected. We measure reflectivities 
as high as 30%. In contrast to a previous experiment devoted to the thin grating limit [S. Slama, et 
al., Phys. Rev. Lett. 94, 193901 (2005)] we now investigate the thick grating limit characterized by 
multiple reflections of the light beam between the atomic layers. In principle multiple reflections give 
rise to a photonic stop band, which manifests itself in the Bragg diffraction spectra as asymmetries 
and minima due to destructive interference between different reflection paths. We show that close 
to resonance however disorder favors diffuse scattering, hinders coherent multiple scattering and 
impedes the characteristic suppression of spontaneous emission inside a photonic band gap. 

PACS numbers: 42.50.Vk, 42.55.-f, 42.60.Lh, 34.50.-s 



I. INTRODUCTION 

The idea of realizing photonic band gaps (PBG) in op- 
tical lattices, i.e. periodic arrays of cold atomic clouds 
confined inside standing light waves, has been published 
in 1995 by Deutsch et al. 0]. A first step towards an 
experimental observation of this phenomenon was made 
by showing that lattices of atomic gases can give rise 
to Bragg scattering in the very same way as X-rays are 
scattered in structure analyses of solid crystals or 
molecules 0- This demonstration has been given by G. 
Birkl et al. and M. Weidemiiller et al. 0, || using near- 
resonant optical lattices. 

In resonant lattices the optical trapping potential pro- 
vides an efficient cooling mechanism for the atomic cloud. 
This cooling is important, because it balances heating 
due to inelastic scattering processes, which destroy the 
periodic order and decrease the Bragg scattering effi- 
ciency. Cooling is absent in conservative optical lat- 
tices tuned far from atomic resonances. On the other 
hand, conservative lattice potentials are interesting in 
view of their perspectives to mimic solid state physics. 
For example, Mott insulator phase transitions in de- 
generate atomic quantum gases have been observed 0, 
and fermionic gases confined in optical lattices are ex- 
pected to exhibit novel quantum phases involving high- 
temperature superfluidity. Those phases may consti- 
tute useful toy models for superconductivity in high-T c 
cuprates Q. Bragg diffraction could represent a novel 
and powerful tool for sensitively probing the properties 
of such optical crystals provided the destructive influence 
of resonant probe light absorption is mastered. 

PBGs are today extensively studied in crystals and 
fibers. Dielectric materials offer the possibility of realiz- 
ing complex periodic structures in three dimensions alter- 
nating high and low index of refraction domains. Those 
structures, called photonic crystals, can exhibit ranges of 
frequencies known as photonic band gaps for which the 



propagation of electromagnetic waves is classically for- 
bidden in certain directions Q. Tailoring of the density 
of states for the electromagnetic modes allows for con- 
trolling fundamental atom-radiation interactions in solid 
state environments and even to suppress vacuum fluctu- 
ations. The hallmarks of a PBG are the inhibition of 
spontaneous emission, an effect that has been observed 
with optical cavities Q and the possibility of Anderson 
localization of light by point defects added to the pho- 
tonic band gap material. 

Although impressive progress has been made [ljj m 
fabricating photonic crystals, they suffer from fundamen- 
tal difficulties in guaranteeing the required fidelity over 
long ranges due to fluctuations in position and size of 
the building blocks. This disorder perturbs those prop- 
erties of photonic crystals based on global interference: 
It reduces the Bragg reflectivity, extinguishes the trans- 
mitted light, and ultimately destroys the photonic band 
gap. On the other hand, optical lattices exhibit an intrin- 
sically perfect periodicity. Local disorder introduced by 
thermal fluctuations in the atomic density distribution 
at each lattice site reduces the value of the Debye- Waller 
factor 12], but does not affect the quality of the long- 
range order. 

To observe photonic band gaps with optical lattices, 
one must reach the thick grating regime. However, all 
Bragg scattering experiments on optical lattices have so 
far [J 0, Il3| been performed in the thin grating regime, 
where the lattice's optical density is so low that multi- 
ple light scattering events are rare. Bragg scattering at 
thin lattices is understood as resulting from constructive 
interference of the Rayleigh-scattered radiation pattern 
emitted by periodically arranged point-like sources. In 
this regime, the reflection coefficient of the lattice turns 
out to be nearly real, phase shifts are negligible, and 
spectral lineshapes are symmetric. In particular, since 
the scattering takes place as a local process, which means 
that the light is scattered by individual atoms, the atomic 
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positions do not shape the absorption spectrum |14| . 

In contrast, the thick grating regime is characterized 
by multiple reflections of the incident light between the 
stacked atomic layers (Bragg planes). The interference 
between the light reflected from or transmitted through 
the layers gives rise to stopping bands for certain light 
frequencies or irradiation angles. In this regime, absorp- 
tion can generally be neglected, but large phase shifts oc- 
cur, and the lineshapes are asymmetric. Multiple beam 
interference globalizes the scattering process. 

In this work we study a one-dimensional optical lattice 
consisting of a standing light wave filled with trapped 
atoms. The atoms arrange themselves into a linear array 
of lens-shaped clouds aligned along their symmetry axis. 
The clouds have a finite radial extent and are centered at 
the locations of the antinodes. We show that for this con- 
figuration the thick grating regime is within experimen- 
tal reach. In fact, which regime is realized in experiment 
depends on the available effective number of scattering 
layers. The number of populated antinodes sets an upper 
limit. However the finite radial extent of the layers also 
limits the number of multiple reflections when the angle 
of incidence for the Bragg light is large. Therefore, to 
get into the thick grating regime, we have modified the 
setup of a previous experiment [13| in order to reduce the 
Bragg angle and increase the radial extent of the Bragg 
layers. With this setup we encounter a new limitation: 
When the probe laser is tuned close to an atomic reso- 
nance, detrimental absorption due to disordered atoms 
reduces the effective number of layers involved in multi- 
ple scattering and dominates the spectra. Off-resonance, 
in contrast, the weak Bragg scattering efficiency brings 
us into the thin grating regime. Hence, it is important to 
identify clear signatures of multiple scattering in optical 
lattices and to develop sensitive tools for their detection. 

In some respects, we may understand the ID optical 
lattice as a dielectric mirror with layers made of a dilute 
atomic gas. On the other hand, as shown in Ref. [l5j . 
the ID optical lattice shares the peculiarities of a linear 
array of point-like scatterers. 

Several authors |4 Ul2tll6| suggested optical lattices for 
the creation of 31? photonic band gaps. However, a con- 
sequence of the narrow width of the atomic resonance is 
that Bragg scattering at gaseous lattices is intrinsically 
one- dimensional. This holds also for 2D and 31? geome- 
tries of optical lattices, so that the results of our investi- 
gations apply to all kinds of lattice configurations. The 
low dimensionality of the scattering problem thus com- 
promises the realization of true 3-D photonic band gaps 
with optical lattices. 

We organized this paper as follows: In section ITT1 we 
present our experimental setup, show examples of typi- 
cal spectra, and discuss how diffuse scattering interferes 
with multiple reflections. To qualitatively understand 
the spectra a transfer matrix model is developed in sec- 
tion |ini Based on Ref. Q the model is extented to com- 
ply with partially disordered lattices and inhomogeneous 



Stark shifts. The model provides a simple picture for 
the observed spectra allowing for a discrimination be- 
tween diffuse and multiple scattering. It also describes 
the expected suppression of spontaneous emission and 
yields a quantitative prediction of the lattice's reflection 
and transmission as a function of experimental parame- 
ters. Section llVI presents our observations and discusses 
them in terms of multiple reflections and diffuse scatter- 
ing. Since both effects give rise to similar signatures in 
the reflection spectra, they are difficult to separate, in 
particular in the presence of experimental imperfections. 
As summarized in section [V] despite the fact that some 
signatures strongly suggest the concurrence of multiple 
reflections, it seems actually beyond reach to see suppres- 
sion of spontaneous emission due to a reduced density of 
optical states. 

II. EXPERIMENT 
A. Setup 

The optical layout of our experiment shown in Figgis 
identical to the one presented in Ref. 0]. It consists of 
an optical cavity and a setup for Bragg scattering. The 
light of a titanium-sapphire laser operating at Adi P = 808 
- 812 nm, which is red-detuned with respect to the rubid- 
ium D\ line, is coupled and phase-locked to the cavity. 
The standing wave which builds up inside the cavity has 
a periodicity of ^Adip = 7r/fcdip- The beam diameter 
at the center of the cavity is u>di P = 220 /im. The in- 
tracavity light power is P cav = 5 W. Between N = 1 5 
and 10 7 85 Rb atoms are loaded from a standard magneto- 
optical trap (MOT) into the standing wave. About 10000 
antinodes are filled with atoms. Typically the tempera- 
ture of the cloud is on the order of few 100 /iK. In ear- 
lier experiments |l5j we found that the temperature of 
the cloud tends to adopt a fixed ratio with the depth 
of the dipole trap, k B T w 0AU o [T^l- Therefore, the 
spatial distribution of the atoms does not vary much 
with the potential depth. From this we derive the rms- 
size of an individual atomic cloud along the cavity axis, 
2(7 z = ^Adi P -\/^BT/2J7o ~ 115 nm in the harmonic ap- 
proximation of the trapping potential. The radial size 
is 2cr r = Wdi P \/kBT/Uo « 140 /im and the mean den- 
sity lies between n = 3x 10 9 and 3 x 10 11 cm~ 3 . For 
the present setup we estimate a Debye- Waller factor of 
/DW = e^.p ff ^0.2 0. 

The light used to probe the Bragg resonance is gen- 
erated with a near infrared laser diode operating at 
Abrg = 780 nm. The laser light is passed through an 
acousto-optic modulator (AOM) and an optical fiber and 
collimatcd to a beam waist of u>brg = 800 /im before 
crossing the dipole trap standing wave under an angle 
of f3i ~ arccos(Abrg/Adi P ), which at Adi P = 810 nm is 
about 15.6°. The incident laser power, P ; — 30 /xW, 
is well below saturation. Some time after loading the 
atoms into the standing wave the probe beam is switched 
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is more than two orders of magnitude higher than in any 
previous measurement on optical lattices. 



FIG. 1: The experimental setup consists of a cavity pumped 
with a titanium-sapphire laser at 808 - 812 nm and a diode 
laser at 780 nm for Bragg scattering. About 10000 antinodes 
of the cavity mode are loaded with 85 Rb atoms from a MOT. 
Observed are the reflected (PMT1), the transmitted (PD) and 
the absorbed (PMT2) light powers. 



on and frequency-ramped across the rubidium D2 reso- 
nance. The light power reflected from the atoms P r is de- 
tected under the angle f3 B = with a photomultiplicr 
(PMT1). The transmitted light power P t is recorded 
with a photodiode (PD), and the isotropically scattered 
power P a is detected by collecting the light emission into 
a solid angle of 0.05 sr orthogonal to the incident probe 
beam with a second photomultiplier (PMT2). 

To obtain Bragg reflection the angle of incidence of the 
probe laser has to be matched to the lattice constant. 
Experimental fine tuning of the Bragg condition is how- 
ever easier by varying the wavelength of the lattice laser, 
AAdip = Adip — Abrg/ cos/3i, while the angle of incidence 
is kept fixed. 



B. Bragg spectra 

The experimentally accessible quantities are the re- 
flected, transmitted and absorbed light powers. We take 
simultaneous spectra of these quantities by ramping the 
probe laser frequency across the Bragg resonance. In or- 
der to compare with calculations we are interested in the 
reflection, transmission and absorption coefficients R, T, 
and A. A direct comparison is complicated by the fact 
that the probe beam cross section is larger than the size 
of the atomic cloud, so that only a fraction 77 w 16% 
of the incident power P really overlaps with the atomic 
cloud H3, 



P r = RriPx , 

P t = T77P + (1 - r))Pi 

Pa = AllP . 



(1) 



Therefore the energy conservation requirement, p = P r + 
P + P a , implies R+ T + A = 1. 

Fig- Efa) shows reflection, transmission and absorption 
spectra of the Bragg resonance obtained by ramping the 
detuning of the probe laser Ab rg from the D2 resonance. 
The resonance linewidth is r/27r = 6 MHz. Note the 
high Bragg reflection efficiency of more than 30%, which 
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FIG. 2: (a) Reflection (lower thick trace), transmission (up- 
per trace) and absorption spectra (noisy thin trace) of Bragg 
scattering, (b) Dependence of the absolute maximum of the 
reflection spectra on the atom number. The maximum Bragg 
diffraction efficiency saturates for large atom numbers. The 
theoretical curve assumes iV s = 400 scattering layers and 
a Debye-Waller factor of /dw = 0.2. (c) Measured Bragg 
spectra for various atomic densities. The small peaks at 
Abrg = — 31P and — 20r corresponds to transitions from the 
ground state hyperfme level F — 3 to the excited state lev- 
els F' — 2, 3, the broad peak around Ab rg = to the level 
F' = 4. For the calculations (d) the optical cross sections at 
the various hyperfine transitions have been weighted accord- 
ing to their relative oscillator strength. 



The scattering efficiency for standard Bragg diffrac- 
tion depends quadratically on the atom number. This 
dependency has been observed in Ref. [13J in the thin 
grating regime. The measurements performed with the 
present apparatus exhibit a different behavior. As seen 
in Fig. Hfb) at large atom densities the scattering ef- 
ficiency seems to saturate. Also the shape of the re- 
flection signal depends critically on the atom number. 
At very low atom number we find a Lorentzian line- 
shape, whose width corresponds to the natural linewidth 
of the D2 line. When the atom number is increased, 
the linewidth broadens and the peak saturates. At large 
atom numbers a pronounced dip appears at the center 
of the strongest resonance of the reflection spectrum Q , 
whose contrast increases with increasing atomic density 
[cf. Fig. GJc)] . Theoretical predictions based on transfer 
matrices calculations presented below confirm this be- 
havior [cf. Fig. Hfd)]. The main goal of this paper is to 
explain these observations. 
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C. Specular versus diffuse scattering 

The atoms are strongly localized in axial direction at 
the center of the antinodes. Consequently, the spectrum 
of the light reflected into the Bragg angle is narrowed by 
the Lamb-Dickc effect, and the axial thermal distribution 
of the atoms does not broaden the angular distribution of 
the reflected radiation, but increases the background of 
isotropically distributed diffuse scattering. This behavior 
is known from Bragg scattering of X-rays at solids. In 
contrast the weak radial confinement of the atoms does 
broaden the angular distribution . 

The finite radial atomic distribution has however a 
more important impact on the reflected light. The ra- 
dial size of the atomic cloud determines the length of 
the probe beam trajectory across the lattice. Thus for 
thick lattices, when multiple scattering plays a role, the 
maximum number of reflections, i.e. the effective num- 
ber of layers, is limited to iV s = 2o>/(Adip tan/?j). In an 
earlier experiment we have studied Bragg scattering 
at an optical lattice near the 5Si/ 2 - 6P3/2 transition at 
420 ran. The corresponding Bragg angle was 58° yield- 
ing N s « 100. In the present setup we use a linear cavity 
with a larger mode waist and operate the probe light near 
the £>2 line. This decreases the Bragg angle to 15.6° and 
increases the effective number of layers to N s ps 600. 

At first glance, the saturation behavior observed in 
Fig. EJb) could be interpreted as an indication for mul- 
tiple scattering. The dips appearing in Fig. |2fc) could 
result from different multiple scattering trajectories of 
the probe laser along the lattice, which destructively in- 
terfere for certain values of the incident angle and of the 
probe light frequency. However, as already demonstrated 
in Ref. [4j the situation is more complex. In fact in a ther- 
mal lattice, for a small Debye- Waller factor, /dw *C 1, 
disordered atoms have a dramatic influence on the scat- 
tering: 1 . They reduce the number of ordered atoms con- 
tributing to Bragg scattering. 2. They absorb and attenu- 
ate the incident light, and hence decrease the penetration 
depth of the probe beam and thus the effective number 
of layers available to Bragg and in particular to multiple 
scattering. 

A first approach in describing the physical situation 
may consist in dividing the atomic cloud into two parts: 
a perfectly ordered optical lattice with density n/ow 
and a homogenous cloud having the density n(l — /dw)- 
The homogeneous cloud limits the penetration depth to 
z p d = [<m(l — /dw)] ) which corresponds to a reduced 
effective number of layers N S!P a = 2z p( j/A. The effective 
number of layers thus critically depends on the density 
of disordered atoms and, via the optical cross section a, 
on the detuning Ab rg - On resonance assuming a den- 
sity n ps 3 x 10 11 cm~ 3 we estimate iV SiP d - 37 0. 
As a consequence, we expect a dramatic break-down of 
the Bragg reflection signal close to resonance, where the 
absorption by the unordered atoms is largest. The fre- 
quency range for absorption, being on the order of the 



natural linewidth T, is much smaller than the width of 
the reflection signal, which explains the appearance of a 
narrow dip in the reflection spectra Fig.|2[c). In particu- 
lar, signatures of multiple scattering will not show up in 
parameter regimes, where the number of scattering layers 
is considerably reduced below N s , N SjP d <C N s . 

Note that the separation of the Bragg reflection into a 
perfectly ordered contribution and diffuse scattering cor- 
responds to an ansatz frequently made in treating disor- 
der and impurities in crystals [12, |2(j ■ 

III. THEORETICAL MODELING 

A. Transfer matrix formalism 

In the following we attempt to provide a more quanti- 
tative understanding of the observations by developing a 
simple theoretical model. We start by relating the atomic 
polarizability a (given in Si-units) to the single-layer re- 
flection coefficient £ via 

C = -nSz — (2) 

r 3 \ Ig 1 

— — TIO Z — 

2 2vr i + 2(A brg - A F )/T ' 

where nSz is the surface density estimated for a homo- 
geneous atomic density n and a layer thickness Sz. To 
account for the presence of several hyperfine transitions 
at frequency detunings A p [cf . Fig. d)] , we build the 
weighted sum of the individual oscillator strengths. 

For the description of the collective influence of the 
atoms on the incident light, we use a generalization of 
the transfer matrix formalism presented in Ref. [l|. The 
generalization concerns two major points: First of all, 
in that reference the probe beam was assumed collinear 
with the optical lattice beams, while in our case the angle 
of incidence is significantly different from 0. In fact the 
deviation of the chosen angle from the Bragg angle consti- 
tutes an additional degree of freedom allowing us to tune 
frequency and quasi-momentum independently. Again, 
in practice we detune the lattice constant rather than 
the angle of incidence. A second generalization consists 
in the inclusion of diffuse scattering into the formalism, 
as detailed in the next section. All theoretical lincshapes 
shown in the figures are obtained from this transfer ma- 
trix model. 

We will now supply the basic ingredients of the model 
skipping the details already reported in Ref. PJ. The 
in- and outgoing field amplitudes of the probe beam at 
any axial location z of the lattice (the model is one- 
dimensional) are labeled E + (z) and E~(z), respectively. 
Their variation from one location to another is described 
by a transfer matrix M, such that (in the complex rep- 
resentation) 
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The procedure consists now in dividing the atomic 
sample into layers. The transfer matrix for interaction 
of the probe light with a single infinitely narrow layer 
of the optical lattice characterized by the surface density 
nSz is 



Ac = 



l + *c «c 

-K 1-iC 



(4) 



The transformation of the field amplitudes between two 
such layers separated by Az is described by 



B Az = 



a ik^ rg Az cos (3j 

o 





-ik\, T gA.z cos /3j 



(5) 



Hence the total transfer matrix for a lattice with N s 
layers reads M = (A^Baz) N " ■ Finally the reflection coef- 
ficient R = 
are calculated via 



and the transmission coefficient T = \t\ 



M 

Mi 



12 



M 2 



(6) 



B. Sequential density model 

Instead of separating the atomic cloud into a perfectly 
ordered lattice and a homogeneous density distribution 
as proposed in Sec. Ill CI we may subdivide every layer 
into a number N ss of sublayers for which we evaluate the 
transfer matrices based on the local density, 



n\oc{z) 



I — 

2N ss ja- u i*')/ k * T dz> 



(7) 



where U(z) = — Uq cos k^ipZ is the trapping potential, 
or in the harmonic approximation U(z) — Uok^ ip z 2 — 
kftTz 2 /2a 2 .. Now we recalculate the local single-layer 
reflection coefficient Q\ oc as in Eq. J2J) and set up the 
total transfer matrix, but using the local density n\ oc (z) 
instead of a homogeneous density n and setting the layer 
thickness to Sz = Adi P /2-/V~ ss , 



M = 



n^(^) b Sz 



(8) 



while the phase shift in reflection follows from <f> = 
arctan(Imr/Rer). If we identify a layer with an antin- 
ode of the standing wave, Az — 5z — Adi p /2, we obtain 
from the Eqs. I© Bragg spectra for the case of a perfect 
lattice, such as those shown in Figs. |3fa-c). 




FIG. 3: The figures (a), (b), and (c) show calculated Bragg 
reflection, transmission and absorption spectra for lattices 
with various numbers of layers N 3 . Here we assumed /dw = 1 
and n = 3 x 10 11 cm -3 . The lattice is slightly detuned from 
the Bragg condition, AAdi p = 0.8 nm, so that the lineshapes 
become asymmetric, (d) Sketch of the photonic dispersion 
relation. The opening of a band gap for values of the quasi- 
momentum 9 near the edges of the Brillouin zone indicates 
a suppression of the density of states p(u>) as drawn in fig- 
ure (e). 




FIG. 4: Calculated Bragg reflection spectra for various detun- 
ings of the lattice constant from the Bragg condition (a) with- 
out and (b) with Stark-shift. We assumed n = 3 x 10 11 cm -3 , 
iV s = 200 and iV ss = 20. At a potential depth U = 500 yiK 
we assume &bT w QAUq. 



At finite temperatures the atoms are distributed over 
the optical potential and thus experience individual dy- 
namical Stark shifts of their resonances varying with 
the atoms' locations. This leads to serious inhomo- 
geneous broadening of the Bragg spectra as shown in 
Ref. 0| . Cold atoms which concentrate at the antin- 
odes of the standing wave potential are Stark-shifted by 
a large amount and form the blue edge of the line pro- 
file. The hot part of the cloud sees on average a shal- 
lower potential and forms the red tail of the profile. A 
possible approach to describe the line broadening con- 
sists in building the convolution of the spectra calcu- 
lated from Eqs. © and © with the probability density 
of finding an atom at a given potential energy. How- 
ever this approach does not account for the fact that in 
the thick grating limit cold (predominantly ordered) and 
hot (mostly disordered) atoms yield qualitatively differ- 
ent lineshapes of the Bragg reflection signal. In other 
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words, the convolution procedure is incompatible with 
the fact that the contributions of cold and hot atoms 
to the Bragg-scattered light depends on the penetration 
depth, which itself varies with the frequency detuning of 
the probe beam. 

Fortunately, the local Stark shift is easily included in 
the sequential densities model via the substitution 



A 



brg 



Abrg - 17(20 



(9) 



The figures E^a,b) show calculated reflection spectra 
without and with Stark shifts for various detunings of 
the lattice constant from the Bragg condition. As in 
Figs. Etb,c) the dip in the spectra corresponds to a joint 
impact of diffuse and multiple scattering and will be dis- 
cussed in more detail in Sec. IfV Al Here we just point 
out that Stark broadening obviously induces pronounced 
asymmetries with respect to AAdip. 

We verified that in the limit of a perfect lattice, ob- 
tained for T — > 0, the line broadening disappears. The 
lines are just blue-shifted by an amount Uq/H. In the 
thin grating regime, obtained for small densities n < 
10 11 cm -3 , the convolution approach and the sequential 
densities approach yield identical spectra. 
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FIG. 5: (a) Evolution of the light intensity across the opti- 
cal lattice normalized to the incident probe beam intensity for 
AAdip = (large amplitude oscillations) and AAdip = 200 nm 
(small amplitude oscillations). The last detuning has been 
chosen to be extreme in order to visualize the rapid decay of 
the intensity even for a low number of layers iV s = 10. The 
probe was tuned to resonance, Abrg = 0. The spikes at the 
bottom of the pictures indicate the maxima of the atomic 
density distribution assumed to describe a perfect lattice, 
/dw = 1- (b) Same as (a) but with /dw — 0.03. The low- 
est curve shows the axial atomic density distribution. (c,d) 
Same as (a) and (b) respectively, but now the probe laser has 
been detuned to Ab rg = F 



C. Suppression of absorption 

The hallmark of 3D photonic crystals is the suppres- 
sion of spontaneous emission. The models describing the 
propagation of light inside photonic crystals assign this 
suppression to a reduction of the density of optical modes 
available for spontaneous decay || . In fact the frequency 
dependence of the density of states, sketched in Fig.^e), 
is characterized by forbidden bands. As pointed out in 
Ref. 0, ID optical lattices exhibit similar band gaps. 
This can be seen in Fig. |3fa-c). For large numbers of 
scattering layers, N s > 1000, the lattice gets opaque. The 
transmission vanishes over a large range of frequency de- 
tunings Abrg, while the reflection is close to unity. Even 
more interesting is the feature that the absorption spec- 
trum 1 — R — T splits into two peaks, and the absorption 
vanishes in the center of the band gap. 

At first glance this may seem surprising. Since the 
geometry of our lattice is ID, we would not expect a 
noticeable modification of the density of decay modes. 
This is similar to the inhibition of spontaneous emission 
inside linear optical cavities [9J: An excited atom may 
decay into all the transverse modes leading out of the 
cavity. If the cavity only covers a small solid angle of the 
radiating atom, its spontaneous decay rate will only be 
reduced by a small amount. 

The reduction of the absorption inside a band gap of 
an optical lattice has a different origin. A deeper un- 
derstanding of the system is gained by calculating the 
progression of the probe light intensity along the optical 
lattice using the above transfer matrix formalism under 
various conditions (cf. Fig.[SJ). First of all, we find that 
the standing wave formed by the incident probe beam 
and the Bragg-reflected light adjusts its phase such that 
its nodes coincide with the atomic layers. In that way 
absorption is minimized. If the length of the lattice is 
finite, the contrast of the standing wave is smaller than 
1, i.e. the probe light intensity at the locations of the 
atomic layers does not vanish. Hence a finite absorption 
subsists even, when the lattice is perfect and the Bragg 
condition fulfilled. 

Let us now study the response of the probe standing 
wave to variations of experimental parameters. Fig. [5fa) 
compares the cases of an aligned and a misaligned angle 
of incidence with respect to the Bragg angle. For a mis- 
aligned angle the periodicity of the probe standing wave 
does not coincide with the lattice constant, which results 
in a displacement of the nodes from the atomic layers 
and hence in enhanced absorption. The contrast of the 
probe wave is smaller than for an aligned angle. However, 
as the probe light penetrates deeper into the lattice, the 
displacement gets smaller and the probe wave contrast 
adopts the value of the aligned case. 

The curves shown in Fig. |5Ta) assume a perfectly 
ordered lattice. In the presence of disordered atoms, 
i.e. atoms which are not confined to the locations of the 
probe beam intensity minima, absorption can take place 
anywhere. We thus obtain an additional background of 
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FIG. 6: Reflection (a) and absorption spectra (b), calculated 
for various values of the Debye- Waller factor /dw. The other 
parameters are AAdip = 0, n — 3 x 10 11 cm -3 , and N s = 600. 
(c,d) Same as (a) and (b), but with AAdip = 0.8 nm. 



absorption, which in the extreme case of strong disor- 
der leads to a fast exponential decay according to the 
Lambert-Beer law. The exponential curves in Fig. |Sfb) 
are obtained under the same conditions as in (a), but 
with a finite Debye- Waller factor, /dw = 0.03. Finally 
when the probe beam is tuned off resonance, the absorp- 
tion is smaller and the penetration depth is drastically in- 
creased. This is shown in Figs.[SJc-d), which corresponds 
to (a-b) respectively, but with a probe laser detuning set 
to Abrg = r. 

In summary, absorption occurs for two reasons: Either 
the Bragg angle is mismatched, or the atoms are dis- 
ordered. This is illustrated by calculations of reflection 
and absorption spectra shown in Fig. El For example, the 
spectrum in Fig. Etd) corresponding to /dw = 1 shows 
an absorption peak despite perfect atomic ordering, be- 
cause the Bragg angle is not matched; and the spectra 
in Fig. EJb) corresponding to /dw < 1 show finite ab- 
sorption although the Bragg angle is matched. Reduced 
absorption is only observed for a perfect lattice and a 
matched angle of incidence for the probe beam [see curve 
corresponding to /dw = 1 in Fig.[(Jb)]. 



IV. OBSERVATIONS 

A. Signatures of multiple scattering 

The appearance of signatures for multiple scattering 
depends very much on the available effective numbers 
of layers. As we have seen in Sec. Ill Cl diffuse scattering 
drastically reduces the number of layers from N s — 600 to 
N S)V d = 37, which is clearly insufficient to produce dips 



in the reflection spectrum corresponding to destructive 
interference of different reflection paths. FigJSfa) reveals 
that for our atomic densities several hundred layers are 
necessary to significantly distort the spectrum near the 
photonic band edge. 
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FIG. 7: Measured (a) and calculated (b) Bragg reflection 
spectra for various dipole laser frequencies. Measured (c) 
and calculated (d) absorption. The calculations are based on 
n = 3 x 10 11 cm -3 , N B = 200, and /dw = 0.7 (dotted lines) 
and /dw = 1 (solid lines). 

In order to observe signatures of a band edge, the angle 
of incidence of the probe beam must be chosen in such a 
way that they appear outside the detuning range, where 
diffuse scattering is dominant. In fact already at |At, rg | > 
3T the penetration depth allows for an effective number 
of layers larger than 600. 

Fig. Eta) shows Bragg reflection spectra measured for 
various lattice constants. Fig. Etc) shows simultaneously 
recorded absorption spectra. Figs. E{b,d) represent cor- 
responding calculations. Although the reflection spectra 
are in qualitative agreement, there is a striking difference. 
While for the measured spectra the dip is always close to 
the line center, the symmetry of the calculated spectra 
changes when going from negative to positive AAdip- An 
explanation for this is given in Sec. II V 51 

The amount of absorption measured in the experiment 
does not depend much on Adi p . This is confirmed by the 
simulations [dotted lines in Fig.EJd)]. The reason for this 
is diffuse scattering due to atomic disorder as discussed in 
Sec. IIII CI pll |22| . For comparison we have also plotted 
in Fig. Efd) the spectra calculated for zero temperature 
(solid lines). The effect of absorption reduction expected 
for /dw — 1 is not observed in experiment. 

Fig. |Sfa) shows a set of spectra recorded for various 
atomic densities. Fig. Elb) shows the corresponding cal- 
culations. The spectra exhibit a well-developed splitting 
of the dip structure. This splitting can not be explained 
by diffuse scattering alone, which means that multiple 
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reflections must play a role. 
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FIG. 8: Occurrence of double dips in the reflection spectra. 
Measurements are shown in (a) . The calculations (b) assume 
AA dip = 0.08 nm, n = 3 x 10 11 cm" 3 , N s = 520, U = 0.6I\ 
and T — 90 fjK. 

The number of free parameters is large, therefore it is 
difficult to pin down the precise value of the experimental 
parameters by fitting. However, there is no realistic pa- 
rameter regime for which our theoretical model predicts 
splittings of the dip without the assumption of multiple 
scattering. Hence, we consider the double dip feature of 
Fig. [8] as the first indication for the existence of a ID 
photonic band gap in an optical lattice. 

B. Experimental side effects 

Various experimental deficiencies can make the quanti- 
tative interpretation of the spectra difficult. First of all, 
when the probe laser is detuned from resonance, its light 
will be bent by refraction when it enters the optically 
thick cloud of atoms. On one hand, this slightly modifies 
the angle of incidence and thus its deviation from the 
Bragg angle. Since this deviation depends on the probe 
beam detuning, the parameters AAdip and Abrg are in- 
tertwined in a complicated way. We have measured the 
reflection angle as a function of Ab rg and found variations 
up to 0.1°, which corresponds to AAdip = 0.4 nm. Ad- 
ditionally, the optically thick cloud focuses or defocused 
the incident beam depending on its detuning 23]. 

A second important issue is the impact of the finite 
radial extent of the atomic layers on the reflection angle. 
In Ref. |Tf| we have shown that, although the lattice con- 
sists of a stack of two-dimensional traps with an aspect 
ratio smaller than cr z /a r ea 10~ 3 , with respect to Bragg 
scattering it behaves more like a chain of point-like scat- 
terers, than like a dielectric mirror. Thus the reflection 
angle is not equal to the angle of incidence, but adjusts 
itself in order to fulfill the Bragg condition. This self- 
adjustment of the Bragg condition impedes a controlled 
tuning of AAdip and explains, why the symmetry of the 
traces in Fig. Eta) does not change when the lattice con- 
stant is varied. 

Furthermore, the transfer matrix model is purely one- 
dimensional. It assumes not only radially infinite atomic 
layers, but also a homogeneous density distribution. In 



reality the radial density distribution is rather Gaussian, 
which implies a variation of the penetration depth with 
the distance from the optical axis. The observed reflec- 
tion spectra thus represent an average of reflection spec- 
tra taken for different optical densities. 

Finally, optical pumping between the ground state hy- 
perfine levels and heating due to resonant absorption may 
occur while scanning the probe beam frequency and dis- 
tort the spectra. These effects can however be accounted 
for in a simple rate equation model, which yields very 
good agreement with our observations. 



V. DISCUSSION AND CONCLUSION 

The observations made with our apparatus clearly 
show features beyond Bragg scattering. These are due 
to two effects: Disorder arising from the thermal distri- 
bution of the atoms introduces strong absorption close 
to resonance, which limits the effective number of layers. 
The second effect is multiple scattering between adjacent 
layers. Even though the above effects tend to distort or 
broaden the features, which are characteristic for multi- 
ple scattering, we find unambiguous signatures of mul- 
tiple reflections. At this stage, it is however difficult to 
exactly quantify the number of layers involved in multiple 
scattering. In any case, for our present parameters the 
absorption spectra do not show any significant reduction 
at resonance, so that the qualification of photonic band 
gap seems not adequate. 

An interesting question concerns the signature of 
atomic ordering in the transmitted and the absorbed 
light. In the thin grating regime, one would not expect 
that the atomic positions influence the absorption: the 
behavior of an absorbing atom does not depend on the lo- 
cation of the other atoms. Moreover, unlike for the reflec- 
tion signal, interference plays no role in forward scatter- 
ing nor in diffuse scattering. Hence there is no signature 
of atomic ordering in T and A. The situation completely 
changes in the thick grating regime, where multiple scat- 
tering between subsequent atomic layers leads to inter- 
ference between the light reflected from or transmitted 
through the layers. The globalization of the scattering 
process leads to interatomic correlation: Now it matters 
how the atoms are arranged. Multiple reflections give rise 
to stopping band gaps for certain ranges of light detuning 
or angle of incidence. It might therefore be more unam- 
biguous to look for signatures of photonic band gaps in 
transmission or absorption spectra. 

To conclude we have shown that long-range spatial or- 
der in atomic clouds can have a dramatic influence on the 
scattering of light. We have extented earlier studies on 
Bragg scattering into the regime of thick gratings charac- 
terized by multiple reflections. Although signatures for 
reduced absorption could not be found due to the fatal 
influence of diffuse scattering, this represents a first step 
towards the realization of photonic band structures in 
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optical lattices. 

Differently from photonic crystals or solid state sys- 
tems, the scattering off optical lattices is weak except 
near atomic resonances. Therefore photonic stopping 
bands are expected to be very narrow. This bears the 
advantage that we can tune the optical density over a 
large range. However, the narrow resonance also implies 
that our system is intrinsically one-dimensional. An ex- 
trapolation to three-dimensional systems seems techni- 
cally demanding, first of all because 3D optical lattices 
have low filling factors of typically 0.01. There are how- 
ever examples of lattices with unity filling factor |24| . 
Bose-condensates in the Mott insulator phase may prove 
useful to guarantee a high and regular occupation of the 
lattice sites Q . On the other hand, the sharpness of the 
resonance results in a very narrow tolerance angle for 
the stop band, which will make it difficult to obtain 3D 
PBGs. Coevorden et al. did numerical calculations 
of the band structure of a 3D optical lattice. To obtain 
a 3D PBG around an atomic transition of frequency u>, 



they had to assume an excessively large spontaneous de- 
cay width, T > O.Olw. 

A major advantage of using ultracold atoms would be 
the total absence of diffuse scattering. Other possible 
technical upgrades include the use of standing waves with 
larger beam waists, thus extenting the radial size of the 
atomic layers, and the choice of smaller Bragg angles, 
which could be done by operating the dipole trap a few 
nm red-detuned from the same transition, the probe laser 
is tuned to. Thus, although it seems today difficult to 
compete with photonic crystals in terms of manipulat- 
ing the propagation of light just by choosing a smart 
arrangement of gaseous atoms, there is much room left 
for improvements. 

We acknowledge financial support from the Landes- 
stiftung Baden- Wurttemberg. 
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